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'  IN  AN  INCOMPEESSIBEE  TURBULENT  BOUNDARY  LAYER 
By  Morris  W.  Eubesin 

SUMMARY 


An  analysis  is  performed  to  determine  the  effects  on  the  rates  of 
convective  heat  transfer  produced  ty  variations  of  the  surface  tempera¬ 
ture  with  distance  along  a  flat  plate.  The  analysis  is  confined  to  the 
case  of  a  low-epeed  turb\ilent  boundary  layer  on  a  flat  plate  oriented 
parallel  to  the  free— stream  velocity.  All  the  physical  properties  of 
the  fluid  are  considered  to  be  constant.  The  temperature  distribution 
in  the  boundary  layer  and  the  local  heat— transfer  rate  at  the  surface 
are  obtained  in  the  form  of  definite  integrals,  with  the  integrands  con¬ 
taining  the  prescribed  surface— temperature  variation  or  its  derivative 
with  respect  to  the  distance  along  the  plate.  Numerical  evaluation  of 
these  integrals  permits  application  of  the  results  to  any  type  of  sur¬ 
face  temperatiare  distribution  that  may  occur  physically. 

The  basic  solution  of  this  analysis,  one  applying  to  a  stepwise — 
discontinuous  surface  temperature  in  which  the  temperature  is  constant 
on  either  side  of  the  discontinuity,  is  correlated  with  a  set  of  exist¬ 
ing  experimental  data.  This  work  is  extended  to  show  that  the  existence 
of  a  stepwise— discontinuous  surface  temperature  upstream  of  a  plug— type 
heat  meter  makes  the  results  obtained  from  this  type  of  an  instrument 
subject  to  large  errors. 


For  convenience,  algebraic  equations  are  presented  for  the  surface 
temperature,  the  surface  temperature  gradient  with  respect  to  the  dis¬ 
tance  along  the  plate,  and  the  local  convective  heat— transfer  rate. 

These  equations  are  in  the  form  of  power  series,  with  positive  exponents, 
in  terms  of  the  distance  along  the  plate.  VThen  one  of  the  quantities, 
siu-face  temperature  or  local  heat  transfer,  is  specified  and  the  coef¬ 
ficients  and  exponents  of  the  particular  power  series  are  evalmted, 
the  other  quantities,  together  with  the  total  heat— transfer  rate,  can 
be  determined  directly  by  substituting  these  coefficients  and  exponents 
into  the  respective  power  series . 
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INTRODUCTION 


Several  experiments  and  analyses  have  been  performed  which  indicate 
that  a  variable  sxirface-temperature  distribution  can  produce  either  a 
marked  increase  or  decrease  in  the  local  and  average  convective  heat- 
transfer  rates  to  or  from  a  surface.  In  reference  1  an  approximate  anal¬ 
ysis  is  described  in  which  the  effects  on  the  convective  heat— transfer 
rates  produced  by  a  stepwise-discontinuous  surface  temperature  on  a  flat 
plate  were  determined.  These  effects  were  determined  from  solutions  of 
the  integral  form  of  the  boundary-layer-energy  equations  for  both  lami¬ 
nar  and  turbulent  boundary  layers.  Further^  these  solutions  were 
restricted  to  the  assumption  of  constant  physical  properties  for  the 
fluids.  It  was  found  that  extremely  large  effects  occur  in  the  region 
directly  downstream  of  the  discontinuity  in  the  surface  temperature. 

In  reference  2,  there  is  described  an  analysis  in  which  the  analyt¬ 
ical  solution  of  the  boundary— layer  equations  was  determined  for  the 
compressible,  laminar  boundary  layer  on  a  flat  plate.  The  surface  temr- 
perature  in  this  analysis  was  represented  by  an  arbitrary  polynomial  in 
terms  of  distance  along  the  plate.  It  was  found  that  even  a  continuous 
surface-temperature  variation  could  produce  large  effects  on  the  convec¬ 
tive  heat  transfer  in  the  laminar  boundary  layer. 

In  reference  3^  there  is  described  an  analysis  in  which  the  solu¬ 
tion  of  the  boundary— layer  energy  equation  was  determined  for  laminar 
flow  on  the  surface  of  a  wedge.  When  the  velocity  distribution  in  the 
boundary  layer  was  assumed  to  be  linear  with  the  distance  normal  to  the 
surface  and  all  the  physical  properties  of  the  fluid  were  assumed  to  be 
constant,  a  solution  of  the  energy  equation  was  obtained  analytically 
for  a  stepwise— discontinuous  surface  temperature.  For  the  case  of  a 
flat  plate  this  solution  gives  results  in  agreement  with  the  results  of 
the  solution  for  the  laminar  boundary  layer  in  reference  1,  This  solu¬ 
tion  (reference  3)  for  the  case  of  a  stepwise-discontinuous  surface  temr- 
perature  was  then  extended  to  an  arbitrary  surface  temperature  through 
the  employment  of  an  integral  solution  similar  in  idea  to  that  of 
Duhamel  (reference  4).  Recently  another  analysis  (reference  5)  was  per¬ 
formed  determining  the  heat  transfer  from  a  body  on  which  the  velocity 
at  the  edge  of  the  boundary  layer  and  the  surface  temperature  vary. 

This  analysis  required  the  same  limiting  assumptions  as  in  reference  3 
although  the  mathematical  details  were  somewhat  different.  Again  it  was 
found  from  these  analyses  that  the  surface— temperature  variation  influ¬ 
enced  the  convective  heat-transfer  rates  considerably. 

It  has  also  been  shown  experimentally  that  a  variation  in  the  sur¬ 
face  temperature  with  distance  along  the  surface  produces  large  effects 
on  the  convective  heat  transfer.  In  reference  6,  an  experiment  is 
described  in  which  the  effect  of  unheated  starting  sections  on  the  aver¬ 
age  heat  transfer  was  determined  from  a  cylindrical  probe  having  a 
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constant  surface  temperature  in  the  heated  region  and  with  its  axis  par¬ 
allel  to  the  free-stream  velocity.  It  was  found  that  the  length  of  the- 
unheated  starting  section  does  influence  the  average  rate  of  heat  trans¬ 
fer  on  the  remainder  of  the  prohe. 

In  the  experiment  described  in  reference  7^  it  was  found  that  val¬ 
ues  of  the  local  heat— transfer  rates  measured  on  a  cone  where  severe 
surface— temperature  variations  occurred  deviated  from  predictions  based 
on  theories  in  which  the  surface  temperature  is  assumed  constant.  These 
results  were  obtained  for  both  laminar  and  turbulent  boundary  layers. 

Eecently  an  experiment  was  performed  to  determine  the  effect  of  a 
stepwise-discontinuous  surface  temperature  on  the  local  heat  transfer 
in  the  turbulent  boundary  layer  of  a  flat  plate ^  reference  8.^  These 
data  show  a  marked  effect  on  the  local  heat  transfer^  and  will  be  cor¬ 
related  with  the  results  of  the  analysis  of  this  paper. 

A  quantitative  determination^  either  experimental  or  analytical, 
is  necessary  to  reveal  the  effect  of  a  continuous  variation  of  surface 
temperature  with  distance  on  the  convective  heat  transfer  in  the  turbu¬ 
lent  boundary  layer.  It  is  the  purpose  of  this  paper,  therefore,  to 
make  this  determination  in  an  approximate  analytical  manner  subject  to 
the  limiting  assumptions  of  a  flat  plate,  of  constant  physical  proper¬ 
ties  of  the  fluid,  and  of  no  frictional  dissipation  of  energy  within 
the  boundary  layer. 


SIMBOLS 


a^^  coefficients  of  power  series,  ^  per  foot"^ 

39+g8m 

A  symbol  defined  by  equation  (22),  foot 
A^  coefficients  of  power  series,  dimensionless 

—1  /35 

B  symbol  defined  By  equation  (42),  foot 

C  coefficient  defined  by  equation  (25),  dimensionless 

c^  specific  beat  at  constant  pressure,  Btu  per  slug,  °F 

F  coefficient  in  equation  (37)^  dimensionless 
1 

The  author  wishes  to  acknowledge  his  indebtedness  to  Mr.  Steve  Scesa 
and  the  Department  of  Mechanical  Engineering,  University  of  California, 
Berkeley,  for  their  kind  peimiissian  to  use  hitherto  unpublished  data 
in  this  report. 
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G  coefficient  in  equation  (37)^  dimensionless 

li  local  heat-transfer  coefficient,  Btu  per  second,  square  foot,  ^  ^ 

h*  local  heat— transfer  coefficient  defined  by  eqmtion  (66),  Btu  per 

second,  square  foot,  ^ 

h  average  heat— transfer  coefficient,  Btu  per  second,  square  foot, 

H  coefficient  defined  by  equation  (40),  dimensionless 
k  summation  index,  dimensionless 

L  distance  along  flat  plate  to  point  of  the  stepwise  discontinuity  of 
the  surface  temperature,  feet 

m  parameter  introduced  in  equation  (20),  dimensionless 

n  summation  index,  dimensionless 

Nu  ITusselt  number  (hx/X),  dimensionless 

P  function  defined  by  equation  (3^)^  dimensionless  ^ 

Pr  Prandtl  number  (iiCp/x),  dimensionless  ^ 

q  local  heat— transfer  rate  per  unit  area,  Btu  per  second,  square  foot 

Q  total  heat— transfer  rate  per  unit  width,  Btu  per  second,  foot 

Be  Reynolds  number  (u^x/v),  dimensionless 
s  total  length  of  heated  portion  of  plate,  feet 
t  temperature  of  fluid  in  the  boundary  layer,  ^ 

u  local  velocity  in  the  boundary  layer  parallel  to  plate,  feet  per 

second 

u^  free-stream  velocity,  feet  per  second 

U  upper  limit  of  integral  defined  in  equation  (^8),  feet 

V  local  velocity  in  the  boundary  layer  normal  to  plate,  feet  per  • 

second 

W  distance  from  leading  edge  of  flat  plate  to  rear  of  heated  section  * 

following  a  surface— temperature  discontinuity,  feet 
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X  distance  along  plate  from  leading  edge,  feet 
y  distance  normal  to  plate  surface,  feet 
function  defined  "by  equation  (65) 
z  temperature  ratio  defined  "by  equation  (38 dimensionless 
5  thickness  of  the  flow  boundary  layer,  feet 
A  thickness  of  the  thermal  boundary  layer,  feet 

Q  exact  temperature  distribution  in  the  boundary  layer  for  a  plate 

at  constant  temperature  preceded  by  an  unheated  starting  section, 
dimensionless 

0*  approximate  temperature  distribution  in  the  boundary  layer  for  a 
plate  at  constant  temperature  preceded  by  an  unheated  starting 
section,  dimensionless 

X  thermal  conductivity  of  fluid,  Btu  per  second,  square  foot,  °F  per 
foot 

p  absolute  viscosity  of  fluid,  pound-seconds  per  square  foot 

V  kinematic  viscosity  of  fluid,  square  feet  per  second 

2  4 

p  mass  density  of  fluid,  pound— seconds  per  foot 
T  local  shear  stress,  pounds  per  square  foot 

Subscripts 


e  effective  property  for  turbulent  flow 
o  free— stream  condition 

w  surface  condition 

p  referring  to  surface  of  heat-meter  plug 

s  referring  to  surface  of  siarrounding  material 


« 
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AMLYSIS 

General  Expressions  for  the  Local  Heat— Transfer  Bate  and  • 

the  Temperature  Distrihution  in  the  Boundary  Layer 


The  method  of  analysis  in  this  paper  will  he  to  determine  an  approx¬ 
imate  solution  of  the  boundary- layer  energy  equation  in  integral  form  for 
a  stepwise-discontinuous  surface  temperature^  and  then  to  extend  this 
solution  to  apply  to  an  arbitrary  surface  temperature  by  employing  an 
Integral  solution  similar  in  idea  to  that  of  Duhamel  (reference  k). 

Boundary— layer  temperature  distribution.—  The  energy  equation  for 
the  turbulent  boundary  layer  on  a  flat  plate^  neglecting  viscous  dissi¬ 
pation,  can  be  expressed  as 


pc 


P 


u  M  +  T 
,  5y  / 


(1) 


where  p,  and  Cp  are  assumed  to  be  independent  of  temperature. 

Suppose  an  exact  solution  of  this  equation  exists  for  the  case  of  a 
plate  which  is  unheated  and  at  the  free— stream  temperature  t^  in  the 
region  x  <  L,  and  is  heated  and  at  the  tempera txare  t^  in  the  region 
X  >  L,  For  the  region  x  >  L 

=  0(x,y,L)  (2) 


where 


e(x,  0,L)  =  0 

e(lH-,y,L)  =  1 

©(ZjOOjL)  =  1 


Because  of  the  linearity  in  temperature  t  of  equation  (l),  various 
solutions  may  be  added  to  satisfy  desired  boundary  conditions.  For 
instance,  the  temperature  in  the  boundary  layer  on  a  flat  plate  which 
is  maintained  at  t„  ,  different  from  t*,  for  x  <  L  and  h,  , 
different  from  for  x  >  L  is  given  by 


(3) 
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in  the  region  x  >  L.  Each  term’ of  equation  (3)  satisfies  the  energy- 
equation  and  the  sum  of  the  terms  satisfies  the  desired  boundary  condi¬ 
tions;  thus,  the  temperature  represented  hy  equation  (3)  can  he  con¬ 
sidered  as  the  formal  solution  of  the  problem.  The  idea  underlying  equa¬ 
tion  (3)  may  be  extended  to  a  large  number  of  temperature  discontinuities 
by  . 


t  —  tij>.  + 


wn  +  ^  L^(k-i)  “  J 


where  x  >  L,  .  Equation  (4)  maybe  rewritten  as 


iJ. 

t  =  t^  0(x,y,O)  +  ^  t^,^  0(x,y,Lj5.)  -  e(x,y,L]5;_i) 


Eqmtion  (5)  is  restricted  to  x  =  L^.  If  n  is  allowed  to  increase 
indefinitely,  the  interval  (Ljj— becomes  small  and  equation  (5) 

may  be  rewritten  in  integral  form  as 


t(x,y)  =  to  e(x,y,0)  +  /  t(L,0)  —  dL 


Equation  (6)  consti-tutes  the  general  expression  for  the  temperature  dis¬ 
tribution  in  the  boundary  layer  for  an  arbitrary  surface— ten^jerature  dis¬ 
tribution. 

Local  heat— transfer  rate.—  The  convective  heat— transfer  rate  per 
\inlt  area  from  a  surface,  where  X  =  Xg,  is  gi-ven  by 


q  =  — X 


In  order  to  introduce  a  temperature  potential  in  the  heat-transfer  expres¬ 
sion,  equation  (6)  is  integrated  by  parts  which  yields 


t(x,y)  =  t(x,0)  -f  [to-t(0,0)]  0(x,y,O)  -  / 


9  dL 
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after  use  of  the  houndary  ccmdltlons  of  eqmtion  (2).  The  local  rate 
of  heat  transfer  obtained  hy  differentiating  and  substituting  equa¬ 
tion  (8)  into  equation  (7)  is 


q(x)  =  -X[tQ-t(0,0)] 

8y 


y=0 


8t(L,0)  ^ 


dL 

y=o 


When  the  local  heat— transfer  coefficient  is  defined  by 


h(x,L) 


X 


a9(x,y^l) 


87 


'y=0 


(10) 


eq-uation  (9)  'becomes 


q(z)  =  h(x,0) 


h(x,L)  dL 

8l 


(11) 


Equation  (11)  constitutes  the  general  expression  for  the  local  heat- 
transfer  rate  per  unit  area.  Note  that  h(x,0)  represents  the  heat- 
transfer  coefficient  on  a  plate  having  a  constant  surface  tomperattu'e. 

It  is  interesting  to  note  that  the  local  '  eat— transfer  rate  per  unit 
area  is  comprised  of  the  local  heat— transfer  rate  on  a  plate  of  constant 
temperature,  equal  to  the  leading-edge  temperature  of  the  plate,  plus 
corrective  terms  introduced  by  the  variation  of  the  surface  temperature. 


Approximate  Solution  for  a  Stepwise- 
Discontinuous  Surface  Temperature 


General  solution.—  The  purpose  of  this  section  is  to  obtain  the 
value  of  h(x,L)  to  be  used  in  equation  (ll).  Becaiise  of  the  complex¬ 
ity  of  the  problem,  an  exact  solution  6,  or  indirectly  h(x,L),  cannot 
be  determined  at  present  and  an  approximate  approach  is  indicated.  Thus, 
the  analysis  presented  in  reference  1  will  be  repeated  herein  in  a  more 
general  fashion. 


The  equation  which  will  be  solved  is  the  boxindary— layer  energy 
equation  a 

V  i  r  , 


O 


(12) 
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The  method,  of  solving  this  eqmtion  is  to  transform  it  to  an  ordinary 
differential  equation  in  terms  of  x  as  the  independent  variable 
and  A  as  the  dependent  variable.  The  resulting  differential  equation 
is  then  integrated. 

Suppose  the  velocity  term  u  in  the  right  member  of  equaticm  (12) 
is  given  by 


where 


(13) 

(14) 


These  are  the  well-known  l/j-power— law  relationships  (reference  9).  In 
addition,  suppose 


even  thou^  6  and  A  may  not  be  of  the  same  order  of  magnitude.  When 
the  velocity  and  temperature  expressed  by  eqmtions  (13)  and  ^5)  are 
substituted  into  the  right  member  of  equation  (12),  there  is  obtained 
after  simplification 


when  A  is  less  than  6'.  It  should  be  noted  that  for  the  case  of  a 
constant  plate  temperature  in  air  A  equals  1.26.  When  the  integra¬ 
tion  of  equation  (12)  is  performed  with  these  limits  on  the  thicknesses 
of  the  boundary  layers  it  is  found  that  the  results  of  the  heat-transfer 
coefficient  differ  from  those  of  equation  (16)  by  the  order  of  0.1  per¬ 
cent.  Therefore,  even  though  A  is  greater  than  6  for  some  cases, 
the  results  of  the  following  analysis  are  believed  to  be  valid  over  the 
entire  plate  to  the  order  of  0.1  percent  for  air. 

In  order  to  transform  equation  (16)  to  an  ordinary  differential 
equation,  it  is  necessary  to  express  the  left  member  of  equation  (16)  in 
terms  of  the  local  characteristics  of  the  boundary  layer.  Many  alterna¬ 
tive  expressions  may  be  obtained  in  this  manner .  This  avoids  imposing 
a  prescribed  variation  of  the  heat— transfer  coefficient  h  in  terms 
of  the  distance  along  the  plate  x,  which  would  make  the  problem 


10 
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trivial*  To  determine  which  of  the  altermtive  solutions  is  applicable, 
correlation  of  a  general  solution  depending  on  a  single  parameter  will 
be  made  with  the  data  of  reference  8.  The  value  of  the  parameter  which 
achieves  the  best  correlation  of  the  data  and  the  analysis  will  then  be 
used  in  the  remainder  of  the  report. 

In  order  to  obtain  these  alternative  expressions  for  h  in  terms 
of  local  characteristics,  it  is  necessary  to  restrict  the  analysis,  for 
the  moment,  to  the  case  of  a  constant  surface  temperature,  For  the  case 
of  a  constant  surface  temperature,  the  Colburn  analogy  between  local 
heat-transfer  rate  and  skin  friction  (reference  10)  is 


(Pr)^''®  = (17) 

P^oCp  pV 

The  ri^t  member  of  equation  (17)  is  expressed,  from  1/7— power— law  con¬ 
siderations,  as 

-1—  =  0.0225  ^  (18) 

P“o 

Further,  it  is  shown  in  appendix  A  that  for  the  case  of  a  constant  sur¬ 
face  temperature  the  ratio  of  the  thickness  of  the  thermal  boundary 
layer  to  the  flow  boundary  layer  is  given  by 

^  (19) 

o 

When  equations  (17) »  (l8)>  and  (19)  are  comhined  so  as  to  express  h/x 
in  terms  of  Pr,  5,  and  A,  the  "variahles  which  appear  in  equation  (l6)^ 
there  restilts 

/  ^3/4  (9— San)  _  i+4m 

^  =  0.0225  Pr  48  A  8™  (20) 

The  thermal  boundary- layer  thickness  was  introduced  into  this  equation 
by  raising  both  sides  of  equation  (19)  to  the  mr-th  power.  Equation  (20) 
constitutes  a  general  expression  for  the  local  heat--transfer  coefficient 
in  terms  of  the  local  boundary— layer  characteristics  and  the  Prandtl 
number.  Although  this  equation  was  derived  for  heat  treuisfer  from  a 
plate  at  a  constant  surface  temperature,  it  will  be  assumed  to  apply  to 
the  case  of  a  flat  plate  with  the  region  x  <  L  unheated  and  at  the 
free-stream  temperature.  It  is  interesting  to  note  that  there  are 
limits  on  the  possible  variation  of  the  parameter  m.  For  instance. 
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in  equaticai  (20)  when  m  is  less  than  -l/4  it  is  found  that  the  heat- 
tr^sfer  coefficient  increases  with  an  increasing  thermal  houndary-layer 
thickness.  When  the  parameter  m  is  greater  than  9/28,  it  is  found 
that  the  heat— transfer  coefficient  increases  with  a  decreasing  Prandtl 
number.  These  variations  are  in  contmdiction  with  existing  knowledge 
concerning  ccaivective  heat  transfer.  The  values  of  the  parameter  m 
must,  therefore,  lie  in  the  interval  -l/4  <  m  <  9/28. 

When  equations  (I6)  and  (20)  are  equated,  there  results  after  sim— 
plification 

2am.+4  ii-s-aam 

=  ""  (21) 

where 


7m^-i 

A  =  (0.2025)(0.37) 


^o\ 

yj 


39+28m  _  39+28111 

140  Pp  48 


(22) 


Eqmtion  (21)  is  of  the  Bernoulli  type  which  can  be  integrated  to  yield 

28 

l9a31+273  39+26811  n'si+S^ 


A  = 


280 


--  L 


280 


280 


)] 


(23) 


when  the  boundary  condition  A  =  0  at  x  =  L  is  imposed.  Equation  (23) 
therefore,  gives  the  thermal  boundary— layer  thickness  on  a  plate  which  is 
unheated  and  at  the  free-stream  temperatiire  to  for  x  <  L  arni  is 
heated  and  at  the  temperature  t^  for  x  >  L.  When  equation  (23)  is 
substituted  into  equation  (20)  there  results 


.  C(m) 


Pr 


1/3 


7+g8m 

39+28m 


(24) 


The  coefficient  in  equation  (24)  is  given  by 

3g  28331—1  T'+SSm 

ri  /  \  \ 00428111  /  39+28331  .  ^^39+28331 

C(in)  =  (0.0225)  (0.37)  (0.0778)  (25) 


In  the  range  -l/4  <  m  <  9/28  the  value  of  C(m)  is  essentially  con¬ 
stant  and  equal  to  0.0288.  It  is  apparent  from  eqmtion  (24)  that  when 
m  is  greater  than  -l/4  and  when  x  approaches  L  the  effect  is  to 
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increase  the  local  heat— transfer  coefficient  many  fold.  When  L  equals 
zero  this  eqmtion  degenerates  to  the  customary  heat-transfer  expression 
based  on  the  1/7— power  laws  for  velocity  and  temperature  distribution 
(reference  11), 

Correlation  with  experiment.—  To  determine  which  value  of  m  makes 
equation  (24)  more  nearly  conform  to  the  measurements,  comparison  is  made 
of  this  solution  with  the  data  of  Scesa  in  reference  8.  The  data  of 
reference  8  were  determined  in  a  well— controlled  experiment  on  a  flat 
plate  oriented  parallel  to  the  free-stream  velocity.  Free-stream  speeds 
of  50  to  70  feet  per  second  were  employed.  The  plate  had  an  unheated 
starting  section  which  was  followed  by  a  heated  section  in  which  there 
was  an  additional  temperature  discontinuity.  All  the  basic  data  except 
those  from  heater  elements  which  can  be  considered  guard  heaters  are 
plotted  in  figure  1.  It  is  observed  that  the  data  lie  consistently  above 
the  theoretical  line  for  a  plate  at  a  constant  temperature. 

The  temperatures  in  the  heated  region  of  the  plate  are  defined  as 
ty^  and  tyg  upstream  and  downstream  from  the  temperature  discontinuity, 
and  the  temperature  of  the  unheated  starting  section  is  t^,  the  free- 
stream  temperature.  The  length  of  the  unheated  starting  section  is  Li, 
and  the  distance  from  the  leading  edge  of  the  plate  to  the  temperature 
discontinuity  in  the  heated  portion  of  the  plate  is  L2.  The  local  heat- 
transfer  rates  per  unit  area  are  then  given,  respectively,  by 


q  =  li(x,Li)(twi-to) 


(26) 


for  Li  <  X  <  La,  and 

q  =  h(x,Li)(t^^-to)  +  li(x,L2)(t^,^-t^,^) 


(27) 


for  X  >  La.  The  heat— transfer  coefficients  in  this  experiment  were 
defined  in  terms  of  the  local  temperature  difference  between  the  surface 
and  the  free  stream.  The  Wusselt  numbers  corresponding  to  those  meas¬ 
ured  in  each  of  these  regions  are  given  by  equation  (24)  as 


eenH-ae 


7+28  m 
39+28m 


Wu  =  (Nu), 


(28) 
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Nu  =  (Nu) 


(i) 


28ni+39  “ 


T+gsm 

39+28  m 


(—') 

\x  / 


ssnH-as  ~ 


7+28m 
39+28  HI 


Jf  the  general  solution,  equation  (24),  yields  results  wMch  conform  to 
experiment,  then  the  experimental  Wusselt  numbers  divided  hy  the  brack¬ 
eted  quantities  should  yield  experimental  Nusselt  numbers  vhich  corre¬ 


spond  to 


These  adjusted  data  are  shown  for  two  values  of  the 


parameter  m,  m  =  0,  and  m  =  9/^8,  in  figures  2  and  3^  respectively. 
It  should  be  noted  that  m  =  — 1/4  corresponds  to  the  data  of  figure  1, 
that  is,  the  bracketed  quantities  are  equal  to  mity.  The  value  m  =  0 
apparently  correlates  the  data  very  well.  This  corresponds  to  the 
results  of  reference  1,  that  is,  the  local  heat— transfer  coefficient  is 
dependent  on  the  thermal  boundary— layer  thickness  eind  the  Prandtl  number 
and  is  independent  of  the  flow  boundary— layer  thickness.  For  m  =  0 
equation  (24)  becomes 

— T/39 

/  ^  1/3  r  /  N,39/40n 

hfeiik  =  0.0288  Pr  J  (30) 


Equation  (30)  will  be  used  in  the  remainder  of  the  analysis. 

Average  heat— transfer  coefficient  determined  from  particular 
solution.—  In  many  cases  it  is  desirable  to  know  the  average  heat- 
transfer  coefficient  over  a  region  at  constant  temperature  preceded  by 
an  unheated  starting  section.  The  total  heat  transferred  in  the  region 
per  unit  width  is  given  by 


q(x)dx  =  h(x,L)(t^-tQ)  dx 


The  symbol  ¥  represents  the  distance  from  the  leading  edge  of  the 
plate  to  the  rear  of  the  region  considered.  The  average  heat— transfer 
coefficient  over  this  region  is  defined  as 


(t^-to)(W-L) 


(32) 
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When  equations  (30),  (31),  and  (32)  are  combined  there  results 


where 


h(W,L)  =  h(W,0) 


4/5 


39  i-(lA^) 


(33) 


(34) 


Equation  (34)  was  evaluated  numerically  and  the  results  are  shown  in 
figure  4, 

The  experimental  data  obtained  by  Jakob  and  Dow  (reference  6)  on 
the  average  heat  transferred  from  a  probe  in  axial  air  flow  preceded  by 
unheated  starting  sections  corresponds,  in  form,  to  the  results  obtained 
from  equations  (33)  and  (34).  It  cannot  be  assumed  tacitly,  however, 
that  the  data  obtained  from  probes  should  compare  identically  with 
theory  based  on  a  flat  plate.  In  particular,  the  data  obtained  on  the 
probe  for  a  constant  surface  temperature  were  about  13  percent  lower 
than  corresponding  theoretical  or  experimental  results  on  a  flat  plate. 
As  any  theoretical  comparison  of  probes  and  flat  plates  indicates  that 
the  rate  of  heat  transfer  from  the  probe  should  be  the  greater,  there 
appears  to  be  some  unexplained  reason  for  the  low  values  of  the  data. 

In  view  of  this,  no  extensive  comparison  of  the  data  and  the  results  of 
the  present  analysis  are  made.  A  single^comparison  for  the  case  of  an 
unheated  starting  section  reveals  that  h/h(W,0)  =  1.37  from  the  Jakob 
and  Dow  results,  while  the  corresponding  result  from  the  present  analy¬ 
sis  is  h/h(W,0)  =  1,45.  This  comparison  is  favorable  when  considera¬ 
tion  is  made  of  the  questions  expressed  previously  concerning  the  com¬ 
parison  of  results  from  probes  and  flat  plates. 


Particular  solution  applied  to  general  single  surface— temperature 
disc ont inui t y ,  —  When  the  region  of  the  plate  preceding  the  tSnperature 
discontinuity  is  at  a  temperature  t^^  other  than  the  free>~stream  temr- 

perature,  the  local  and  average  heat— transfer  coefficients  in  the  region 
L  <  X  <  W  at  a  temperature  t^^  are  given  by 


h(x,L)  =  h(x,0) 


.  +  '  '  . . 

4-4-  4-4- 


-ro' 


39/40- 

—7/39 

/ 

(35) 
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S(W,1)  =  l.(W,0)  5  Lml/V)  J  ^  ,  iiO  a/V).  ^^P(L^) 

^  i-(l/w)  tv2-to  39  i-(l/w)  t^,2-to 


For  certain  applications  in  a  snail  region  (W— L)  it  is  interesting  to 
know  the  ratio  of  the  average  heat— transfer  coefficient  to  the  local 
heat— transfer  coefficient  which  would  premil  at  the  center  of  the  region 
if  the  plate  were  at  a  constant  temperature.  This  ratio  is  obtained  hy 


dividing  equation  (36)  hy  h 


h(W^L)  ^  -p.  ^  .  p  /l\ 

h(^,o)  '  Wt,^-to 


where  F(l/w)  and  G(l/w)  were  evaluated  numerically.  If 


Sis  '^1 


equation  (37)  may  he  rewritten  as 


where 


.  =  F(L/W)  +  H(L/W)z 

-C-fA 


h(l/w)  =  g(l/w)  -  f(l/w) 


The  numerical  values  of  F(l/W)  and  H(l/w)  are  plotted  in  figure  5 


Extension  of  Approximate  Solution  to  the  Problem 
of  an  Arbitrary  Surface  Temperature 


Boundary— layer  temperature  distribution.—  The  approximate  temper- 
atirre  distribution  determined  from  equations  (15)  and  (23)  with  m  =  0 
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0*  = 


BZ 


1/7 


( 


312/280  273/280  39/280  \41/39 

X  -L  X  J 


(41) 


where 

B  =  1.15 


(42) 


At  the  outer  edge  of  the  houndary  layer  0*  becomes  equal  to  unity; 
therefore,  equation  (4l)  is  valid  only  in  the  region 


^  0 


312/280  273/280  39/2e0\a8/39 


) 


For  values  of  y  greater  than  given  by  the  inequality  (43) 


(43) 


0‘  =  1 


(44) 


For  the  range  of  y  indicated  by  the  inequality  (43) 


S0» 

Sl 


10 


1/7  39/280  1/40 


312/280  273/280  39/280\  43/39 

X  -  L  ^  J 


For  greater  y 


(45) 


30 » 

3L 


=  0 


(46) 


When  equations  (4l),  (45),  and  (46)  are  substituted  into  equation  (6) 
there  results 


t(x,y)  =  By 


1/7 


39/280  -1/40 

t(L,0)  X  L  '  dL 

4/35  10  J 

(x  o 

/  312/280  _  j^273/2a0  ^39/28o'\  43/39 

(4?) 
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where 


273/280 


g39/4  y 39/28  ^’-39/280 


280/273 


When  By  =  x  it  is  apparent  from  equation  (48)  that  U  =  0,  and 

consequently  from  equation  (47)  t(x,y)  =  tQ.  It  is  not  obvious,  however, 
that  the  left  member  of  equation  (47)  equals  the  prescribed  surface  temr- 
perature  at  y  =  0.  It  is  shown  in  appendix  B  that  equation  (47)  does 
satisfy  the  boundary  conditions  at  y  =  0. 

I 

Local  beat— tra.ns f fir  The  general  expression  for  the  local 

heat— transfer  rate  is  obtained  by  substituting  the  local  heat— transfer 
coefficient  given  by  equation  (30)  into  equation  (ll) 


q  =  h(x,0)  t(0,0)  -  t.  + 


V  39/40n  — T/39 

X  )  “'I 


For  an  arbitrary  surface-temperature  distribution  it  is  necessary 
to  evaluate  the  integrals  in  equations  (47)  and  (49)  numerically.  There 
is,  however,  a  type  of  surface  distribution  for  which  the  local  heat- 
transfer  rate  can  be  determined  easily  by  analytical  treatment.  This 
will  be  shown  in  the  next  section. 

Surface— temperature  distribution  represented  by  nower  series.—  Let 
the  surface  temperature  be  represented  by 


t(x,0)  -  to  =  )  a^x^  n>  0 


where  n  is  not  necessarily  an  integer.  When  eqmtion  (50)  is  differ¬ 
entiated  with  respect  to  x  there  results 


5t(x,0 

5x 


na^x 


The  substitution  of  equation  (51)  at  x  =  L  into  equation  (49)  yields 


q  =  h(z,0)  {  ao  + 


n— 1 

2^  mnh 


\  39/40-1— V39 

L  ]  ■  dL 


This  equation  may  be  rewritten  as 


q  =  h(x,0)  ao  +  )  ha 


h— 1  /  L 


39/40  _7/39 
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,  /  .39/40 

Letting  a  =  (L/x)  allows  transformation  of  the  integral  terms  of 

equation  (53)  so  that  the  equation  reads 


q  =  h(x,0) 


8-0 


39 


m^x 


n 


1  .iSa— 1 


a' 


39 


-7/39 

(l-a)  da 


n 


(54) 


The  integral  terms  of  this  equation  are  Euler ian  integrals  of  the  first 
kind  (reference  12)^  the  values  of  which  are  given  hy  the  Beta  function 
which  may  he  expressed  in  Gamma  functions  as 


where  n  is  greater  than  or  equal  to  zero.  The  values  of  the  Gamma 
functions  may  he  obtained  from  tables  in  reference  I3.  In  view  of  equa¬ 
tion  (55 equation  (5^)  becomes 


When  the  recursion  equation 


r(y+l)  =  yr(y) 


(57) 


is  used,  equation  (56)  may  be  rewritten  as 


It  is  apparent  that  once  the  type  of  surface— temperature  distribution 
is  specified,  that  is,  values  of  specified,  the  distribution  of 
the  local  heat-transfer  rate  per  *unit  area  can  be  determined  directly 
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Ttie  term  is  plotted  as  a  function  of  n  in  figure  6.  If  the  heat- 

transfer  coefficient  for  the  case  of  a  variable  surface  temperature  is 
defined  in  terms  of  the  local  temperature  difference 


* 

h 


q 

t(z,0)  -  to 


(66) 


then 


(67) 


Inspection  of  the  right  member  of  equation  (67)  reveals  several  inter¬ 
esting  facts.  For  a  problem  in  which  the  denominator  has  a  local  value 
of  zero,  the  numerator  will  be  finite  in  all  but  a  very  exceptional  case; 
thus  h*/h(x,0)  will  be  infinite  at  this  point.  For  another  case  the 
niomerator  may  be  zero  locally  while  the  denominator  will  probably  be 
finite;  thus  h*/h(x,0)  will  be  zero  at  this  point.  It  is  seen,  there¬ 
fore,  that  h*/h(x,0)  can  attain  any  numerical  value,  depending  on  the 
problem  considered. 

In  summary,  it  is  noted  from  equations  (61),  (62),  and  (63)  that  the 
values  of  the  coefficient  may  be  determined  by  a  prescribed  distri¬ 

bution  of  surface  temperature,  a  prescribed  distribution  of  the  surface 
temperature  gradient,  or  a  prescribed  distribution  of  the  local  heat- 
transfer  rate  per  unit  area.  Once  the  coefficients  are  evaluated 

from  the  prescribed  qmntity  the  other  two  quantities,  together  with  the 
total  heat-transfer  rate,  can  be  determined  directly. 


EXAMPLES 

Plug— Type  Heat  Meter 


Convective  heat— transfer  rates  at  a  surface  are  often  measured  con¬ 
veniently  by  determining  the  heat— trans fer  rate  to  or  from  a  plug  which 
is  thermally  isolated  from  the  surrounding  material.  The  plug  can  be 
used  as  a  steady— state  device  in  which  heat  is  generated  electrically 
and  the  temperature  of  the  pliig  is  maintained  constant.  The  usual  prac¬ 
tice,  however,  is  to  use  the  plug  as  a  transient  device  by  measuring  the 
temperature— time  history  of  the  plug.  The  local  heat-transfer  rate 
throu^  the  surface  of  the  plug  is  represented  by  the  product  of  the 
rate  of  change  of  ten5)erature  and  the  thermal  capacity  of  the  plug. 
Regardless  of  the  convenience  of  using  such  a  device,  the  present  ■ 
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analysis  Indicates  that  the  measurements  are  subject  to  certain  inherent 
errors  in  addition  to  those  which  may  arise  from  imperfect  thermal  iso¬ 
lation. 

It  is  not  possible  to  design  a  simple  plug— type  heat  meter  which 
will  always  maintain  a  temperature  identical  to  the  temperature  of  the 
surroundii^g  material.  If  an  extremely  thin  insulating  material  is  used 
between  the  plug. and  the  surrounding  material,  a  surface— temperature  dis¬ 
continuity  will  occur  at  the  seam  of  the  plug  and  the  surrounding  mate¬ 
rial.  If  the  instantaneous  surface  temperature  of  the  plug  is  defined 
as  tp,  while  the  surface  temperature  of  the  surrounding  material  to  the 
leading  edge  of  the  plate  is  denoted  tg,  it  is  apparent  that  eqmtions, 
(36)  through  (4o)  are  applicable  to  this  problem,  where  tp  =  t^^ 

“^s  =  Wi* 


It  is  desirable  to  know  how  well  a  plug-type  heat  meter  measures 
the  local  heat— transfer  coefficient  which  would  exist  at  the  position  of 
the  center  of  the  plug  if  the  plate  were  at  a  constant  temperature. 
Equation  (37)  is  directly  applicable  to  this  problem.  For  example,  sup¬ 
pose  the  instantaneous  temperatures  of  a  heat-meter  plug,  the  surrounding 
material,  and  the  free  stream  are  125°  5*^  100°  F^  and  200^  F,  respec¬ 
tively,  The  dimensions  of  the  plug  are  L  =.24  inches  and  W  =  25  inches. 
It  is  found  from  equation  (38)  that  z  =  -0.33>  and  l/W  =  O.96.  From 
equation  (39)  and  figure  5 


h(25,24) 

11(24.5,0) 


0.62 


(68) 


Thus,  the  plug— type  heat  meter  measures  a  heat— transfer  coefficient  which 
is  only  62  percent  of  the  local  heat— transfer  coefficient  which  would 
have  existed  at  the  position  of  the  center  of  the  plug  had  the  plate  been 
at  a  constant  surface  temperature. 

Although  this  example  is  arbitrary,  the  conditions  are  by  no  means 
implausible.  It  should  be  noted  from  equation  (39)  that  more  severe 
surface-temperature  conditions  can  result  in  negative  values  of 

— — ,  which  means*  that  measurement  of  heat  transfer  by  a  plug 
h(W+L/2,0)  ^  ^ 

Would  be  in  a  direction  opposite  to  that  which  would  normally  prevail  on 
a  plate  having  a  constant  surface  temperature.  These  results  indicate 
that  extreme  caution  should  be  exercised  in  the  interpretation  of  the 
data  obtained  by  plug— type  heat  meters. 
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Constant  Heat— Transfer  Hate 


Suppose  it  is  desired  to  maintain  a  surface— temperature  distritu— 
tion  such  that  the  local  heat— transfer  rate  per  unit  area  is  constant 
along  a  plate.  From  equation  (30)  it  can  be  seen  that  h(x',0>  is 
inversely  proportional  to  the  l/5-power  of  x.  Equation  (63)  can  there¬ 
fore  be  written  as 


q  =  t„  h(s,0)  (^1  Y  \(i'j  (69) 


It  is  apparent  that  to  maintain  q  constant,  n  must  have  the  value 
1/5  only.  Eqmtion  (69)  becomes 


q  =  to  h(s,0)  Ai/5  Ti/5 


(70) 


Therefore, 


1/5 


"to  ^2./b 


(71) 


The  required  surface— temperature  variation  as  given  hy  equation  (6l)  is 


t(x,0) 


1  + 


q _ 

to  Vs 


(72) 


Further 


h(x,0) 


1.06 


(73) 


This  latter  equation  indicates, 
local  heat— transfer  coefficient 
transfer  rate  differs  from  that 
hy  only  6  percent. 


from  the  values  of  figure  6,  that  the 
on  a  plate  having  a  constant  heat- 
on  a  plate  at  a  constant  temperature 
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Surface  Temperature  Eepresented  hy  Polynomial 
Suppose  the  temperature  distribution  on  a  flat  plate  is  given  by 

.  1.5  +  0.5(1)  '  (74) 

Then  in  equation  (61),  n  =  0  and  n  =  1,  and  =  0.5  and  Ai  =  0.5. 

From  equation  (63)  and  figure  6,  the  local  heat— transfer  rates  are  giTOn 
tiy 


- ^ —  =  0.5  +  0.612 

to  h(x,0) 


(75) 


The  ratio  of  the  local  heat— transfer  coefficient  to  the  heat— transfer 
coefficient  which  would  exist  if  the  plate  were  at  a  constant  surface 
temperature  is 


h*  ^  0.5  +  0.612(x/s) 
h(x,0)  0.5  +  0.5(x/s) 


(76) 


The  values  of  the  members  of  equations  (7^)  and  (75)  are  plotted  in 
figure  7*  It  can  be  observed  that  althou^  the  temperature  variation 
along  the  plate  increased  the  over-all  temperature  difference  between 
the  plate  and  the  free  stream  by  100  percent^  the  local  heat— transfer 
coefficient  differs  from  that  on  a  constant— temperature  plate  by  a  maxi¬ 
mum  of  about  11  percent. 


CONCLUSIONS 


The  conclusions  of  this  report  are  subject  to  certain  limitations 
inherent  in  the  assumptions  of  the  aiialysis.  The  following  assumptions 
have  been  used: 

1.  All  physical  properties  of  the  fluid  are  constant,  . 

2.  Frictional  dissipation  of  energy  within  the  boundary  layer  is 
negligible. 

3-  The  velocity  distribution  in  the  boundary  layer  is  of  the 
1/7— power  form. 


2k 
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The  temperature  dis  trihut  ion  in  the  houndary  layer  is  of  the 
f“orm  even  when  the  thermal  boundary— layer  thickness  differs 
considerably  from  the  flow  boundary— layer  thickness. 

5.  The  local  heat— transfer  coefficient  determined  on  a  plate  having 
a  constant  surface  temperat\ire  applies  to  a  plate  having  a  variable  sur¬ 
face  temperature  when  it  is  expressed  by  an  eqmtion  based  on  the  local 
flow  and  thermal  boundary-layer  thicknesses  instead  of  the  distance  along 
the  plate. 

Although  the  analysis  based  on  these  assumptions  is  correlated  with 
experimental  data^  this  is  not  a  verification  of  the  individual  assump¬ 
tions  of  the  analysis,  but  rather,  is  a  justification  of  the  use  of  the 
end  results. 

The  foregoing  analysis,  together  with  the  examples  cited,  has  indi¬ 
cated  that  a  variation  in  the  surface  temperature  with  distance  along  a 
flat  plate  influences  the  local  convective  heat  transfer  to  an  extent 
which  depends  on  the  type  of  variation.  In  general,  a  sudden  surface- 
temperatiire  j-ump,  or  discontinuity,  produces  extremely  large  increases 
or  decreases  in  the  convective  heat  transfer  directly  downstream  of  the 
position  of  the  discontinuity.  From  this  it  can  be  concluded  that  meas¬ 
urements  by  a  heat  meter  of  the  plug  ti^pe  can,  in  most  instances,  deviate 
considerably  from  the  heat  transfer  which  would  normally  exist  at  the 
location  of  the  instrument.  A  continuously  variable  surface  temperature 
was  shown,  in  general,  to  have  a  smaller  effect  on  the  convective  heat 
transfer  than  a  sudden  temperature  discontinuity.  For  instance,  it  was 
shown  that  the  local  heat— transfer  coefficient  on  a  plate  having  a  con¬ 
stant  heat— transfer  rate  along  the  surface,  with  the  necessary  l/5— power 
variation  of  the  surface  temperature,  differs  from  that  on  a  plate  with 
a  constant  surface  temperature  by  only  6  percent.  For  the  case  of  con¬ 
tinuously  varying  surface  temperatures,  large  effects  on  the  local  heat- 
transfer  coefficient  are  expected  only  where  the  surface  temperature 
approaches  that  of  the  free  stream  at  a  point  other  than  the  leading 
edge  of  the  plate. 
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APPEaroix  A 

RATIO  OF  BOUKDART-IAIER  THICK2IESSES  OR  PLATE 
WITH  CONSTAHT  SURFACE  TEMPERATURE 


If  tlie  1/7-power— law  velocity  distribution 


(Al) 


is  substituted  into  the  von  Karm^  momentum  integral  for  a  flat  plate 


there  results 


(A2) 


(A3) 


When  the  indicated  operations  are  performed^  there  results 

=  X  ^  •  (a4) 

puo^  72  dx 


Similarly^  if  the  1/7— power  temperature  distribution 


(A5) 


and  equation  (Al)  are  inserted  into  the  integral  form  of  the  energy 
eq^uation 


u(t-to)  dy 

t^-to  dx  J 


(A6) 


there  results 
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when  is  constant  and  it  is  postulated  that  6  is  greater  than  A. 

For  a  constant  surface  temperature  the  Colhum  modification  of  the 
Reynolds  analogy  is  expressed  as 


h 

PCp  Uq 


2/3 

(Pr) 


T 

pV 


(A8) 


When  equations  (a4)  and  (A7)  are  substituted  into  equation  (a8),  there 
is  obtained 


dx 


(A9) 


Therefore 


^8/7 


s  +  c 


(AlO) 


For  the  case  of  a  constant  surface  temperature,  both  the  thickness  of 
the  thermal  and  flow  boundary  layers  are  zero  at  the  leading  edge  of  the 
plate;  thus,  C  equals  zero  in  equation  (AlO)  and 


A 

6 


(Pr) 


7/lS 


(All) 
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t(x,y)  = 


_  to 


4/35 


M 

39 


^39/4^n+39/28 


>,^40/39 

39/40  —39/280  -  ^ 

X  -X  '  ,0 


^312/280^11 


43/39  —431/39 

cp  y 


When  n  =  —39/28,  equation  (B7)  reduces  to 


t(z,y)  = 


t 

4/35 


39/4 

4b  o"'  t 

I 

312/280 


4b 

39  J 


( 


39/40  —39/280  39/28 

X  -X  y 


>,40/39  - 

cpy  ,o_ 


9 


43/39 


39/28 


At  y  =  0,  equation  (b8)  Becomes 


t(x,0)  =  -  —  P 

39  i 


39/4 


t(x,o)dcp 

43/39 

<P 


On  integration,  equation  (B9)  becomes 


t(x,0)  =  M  t(x,0)  [  ^  q>’ 


/39 


B 


00 


39/4 


which  reduces  to  the  identity 


t(x,0)  =  t(x,0) 


y-^+i/^dcp 

(BT) 


—  (B8) 


(B9) 


(BIO) 


(Bll) 
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Figure  2.  -  Seesaws  data  altered  by  approximate  theory  with  m 


2000 


MCA  TN  2345 


33 


Figure  3,  -  Scesa*s  data  altered  by  appraximate  theory  with  m 
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Figure  5.  -  Numerical  values  of  coefficients  in 


equation  (39). 
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Figure  7.  —  Effect  of  variable  surface  temperature 
the  heat  -  transfer  characteristics  in  example 
problem. 
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